Lyapunov Control on Quantum Open System in Decoherence-free Subspaces 
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A scheme to drive and manipulate a finite-dimensional quantum system in the decoherence- 
free subspaces(DFS) by Lyapunov control is proposed. Control fields are established by Lyapunov 
function. This proposal can drive the open quantum system into the DFS and manipulate it to any 
desired eigenstate of the free Hamiltonian. An example which consists of a four-level system with 
three long-lived states driven by two lasers is presented to exemplify the scheme. We have performed 
numerical simulations for the dynamics of the four-level system, which show that the scheme works 
good. 

PACS numbers: 03.65.-w, 03.67.Pp, 02.30.Yy 



I. INTRODUCTION 

Manipulating the time evolution of a quantum system 
is a major task required for quantum information pro- 
cessing. Several strategies for the control of a quantum 
system have been proposed in the past decade pj, vifhich 
can be divided into coherent and incoherent control, ac- 
cording to how the controls enter the dynamics. Among 
the quantum control strategies, Lyapunov control plays 
an important role in quantum control theory. Several pa- 
pers have be published recently to discuss the application 
of Lyapunov control to quantum systems [2-la|. Although 
the basic mathematical formulism for Lyapunov control 
is well established, many questions remain when one con- 
siders its applications in quantum information process- 
ing, for instance, the Lyapunov control on open system 
and the state manipulation in its decoherence-free sub- 
space. 

As a collection of states that undergo unitary evolu- 
tion in the presence of decoherence, the decoherence-free 
subspaces (DFS) [0] and noiseless subsystem(NS)i7ll are 
promising concept in quantum information processing. 
Experimental realizations of DFS have been achieved 
with photons Q and in nuclear spin systems Q- A 
decoherence-free quantum memory for one qubit has 
been realized experimentally with two trapped ions [lOJ, 
lll| . An in-depth study of quantum stabilization prob- 
lems for DFS and NS of Markovian quantum dynamics 
was presented infl^. 

Most recently, we have proposed a scheme to drive 
an open quantum system into the decoherence-free 
subspaces [g. This scheme works also for closed quantum 
system, by replacing the DFS with a desired subspace. 
The result suggests that it is possible to drive a quantum 
system to a set of states (for example, the DFS in the 
paper), however it is difficult to manipulate the system 
into a definite quantum state in the DFS. The aim of 
this paper is to design a Lyapunov control to drive an 
open system to a definite state in the DFS. The Lya- 
punov control has been proven to be a sufficient simple 
control to be analyzed rigorously, in particular, the con- 
trol can be shown to be highly effective for systems that 



satisfy certain sufficient conditions, which roughly speak- 
ing are equivalent to the controllability of the linearized 
system. In Lyapunov control, Lyapunov functions which 
were originally used in feedback control to analyze the 
stability of the control system, have formed the basis for 
new control design. By properly choosing the Lyapunov 
function, our analysis and numerical simulations show 
that the control scheme works good. 

This paper is organized as follows. In Sec. |lll 
we present a general analysis of Lyapunov control for 
open quantum systems, Lyapunov functions and control 
fields are given and discussed. To illustrate the gen- 
eral formulism, we exemplify a four-level system with 2- 
dimensional DFS in Sec. IIIIl showing that the system 
can be controlled to a desired state in the DFS by Lya- 
punov control. Finally, we conclude our results in Sec. 
lYl 



II. GENERAL FORMULISM 

We can model a controlled quantum system either by 
a closed system, or by an open system governed by a 
master equation. In this paper, we restrict our discussion 
to a TV-dimensional open quantum system, and consider 
its dynamics as Markovian and therefore the dynamics 
obeys the Markovian master equation (?i = 1, throughout 
this paper). 



p = -t[H,p] + C{p), 



(1) 



-[L,nP,L}J), H = 
., M) are positive 



where C{p) = \ Em=i \n{[L„^, pL]^] 
Ho + ELi Ut)H,,. A„,(m = 1,2,, 
and time-independent parameters, which characterize the 
decoherence. L„i(rn = 1,2,...,M) are jump operators. 
Hq is a free Hamiltonian and -ff„(n = 1, 2, ..., F) are con- 
trol Hamiltonian, while /„(t)(n — 1,2,...,F) are control 
fields. Equation ([I} is of Lindblad form, this means that 
the solution to Eq. ([!]) has all the required properties of 
a physical density matrix at all times. 

By definition, DFS is composed of states that undergo 
unitary evolution. Considering the fact that there are 
many ways for a quantum system to evolve unitarily, 



we focus on the DFS here that the dissipative part C{p) 
of the master equation is zero, leading to the following 
conditions for DFS[l^. A space spanned by %dfs = 
(I'/'i), IV'2), •••, IV'd)} is a decoherence-free subspace for 
all time t if and only if (1) T-Ldfs is invariant under 
Hq; (2) Lm\ilM) = Cm\ipn) and (3) T\ipn) = g\ipn) for aU 
n — 1,2, ...,D and m — 1,2, ...,M with g = J2i=i -^/Ic/P, 
and r — J2m=i ^niLl^Lfn- With these notations, the goal 
of this paper can be formulated as follows. We wish to ap- 
ply a specified set of control fields {fi{t), n = 1,2, ..., F} 
in Eq. ([T]) such that p(t) evolves into a desired state in the 
DFS and stays there forever. In contrast to the conven- 
tional control problem J^], we here develop the control 
strategy to open system. 
We use 



V{p) = Tr(pi) 



(2) 



as a Lyapunov function, where A is hermitian and time- 
independent. First, we analyze the structure of critical 
points for V{p) with restriction Tr(p) = 1. To determine 
the structure of V{p) around one of its critical points, 
for example pc = '^jP'j\A-j){Aj\ , wc consider a finite 
variation Sp such that Tr(pc + Sp) — 1- Here wc denote 
the normalized eigenvectors and eigenvalues of v4 by \Ai) 
and Ai {i = 1, 2, 3, ..., N), respectively. Express {pc + Sp) 
in the basis of the eigenvectors of A, 



Pc 



^P = IIpjI^J + 5Aj){Ai + SAjl 



N 



(3) 



Q = l 

The normalization condition Tr(pc + 5p) = 1 follows, 



Then 

v{p, + 5p) ~ v{p,) = Y.'p", E(^" - ^i)^':^'o.- (4) 

Considering (5^ as variation parameters and noting 
S^a^a > 0, we find that the structure of V{p) around the 
critical point pc depends on the ordering of the eigen- 
values: Pc is a local maximum as a function of the vari- 
ations (5^ if and only if Aj is the largest eigenvalue, a 
local minimum iff Aj is the smallest eigenvalue and a 
saddle point otherwise. This observation leads us to sus- 
pect that the minimum of V is asymptotically attractive, 
in other words, the control field based on this Lyapunov 
function would drive the open system to the eigenstate 
of A with the smallest eigenvalue. We will show through 
an example that this is exactly the case. 

Now we establish the control fields f,i{t). V{p) — 
Tr(pi) yields, 

V = Tr{C{p)A) - tTr{p[A, ^ /„(f)i/„]). 



where we choose [A, Hq] = 0, because (a, b any opera- 
tors) Tr[a, b] = 0, i.e., the commutator can never be sign 
definite. The choice of [A, Hq] — implies that Hq and A 
must have the same eigenvectors, then the control field 
would drive the open system into an eigenstate of the 
Hamiltonian Hq. To make T^ < 0, we choose a fjg{t) 
such that 



,^, _ . Tr(/:(p)i) 
^"^^" Tr([i,if,Jp)' 

f,it) = -t^,iTT{[A,H,]p)y 



for j 7^ jo- (5) 



Here kj > will be refereed as the strength of the con- 
trol. Then the evolution of the open system with Lya- 
punov control can be described by the following nonlinear 
equations 



m 



-i[Ho + J2Ut)Hn,p{t)]+C{p), (6) 



where /„ (i) is determined by Eq. ^ . It should be empha- 
sized that fjg always exists. To find fjg, TT{[A,Hjg]p) ^ 
is required. This can be done by construction. Now 
we show that fjg is real. By the definition of C{p), C{p) 
is hermite, then Tr(£(p)yl) can be treated as the time 
derivative of {A) and thus is real. Identifying A with a 
hermitian operator for a system described by the Hamil- 
tonian Hjg, we have i^ = [A,Hjg], so TY{i[A,Hjg]p) is 
real. By the same virtue, we can show that all the con- 
trol fields are real as long as the control Hamiltonian Hj 
(j = 1,2,3...) are hermitian. 

By the LaSalle's invariant principle .IS], the au- 
tonomous dynamical system Eq.® converges to an in- 
variant set defined by 5 = {V = 0}. This set is in general 
not empty and of finite dimension, indicating that it is 
easy to manipulate an open system to a set of states 
but diflacult to control it from an arbitrary initial state 
to a given target state. Fortunately, by elaborately de- 
signing the control Hamiltonian and the operator A, we 
can solve this problem as follows. The invariant set de- 
fined by 5 = {y = 0} is an intersection of all sets £j 
{j = 1, 2,3, ...), each one satisfies, 

TriAH^p - HjAp) = 0, 

leading to [A, p] = 0, [A,Hj\ = or [Hj,p] = 0. By 
elaborately choosing Hj {j = 1,2,3,...), we can set the 
contribution of [^,i?j] = and [Hj,p] = to the inter- 
section (i.e., the invariant set £) to zero. In this case, 
the invariant set is a collection of state {pin} that satis- 
fies [A, Pin] = 0. Considering that only the states in DFS 
are stable, we claim that we can manipulate the system 
from any initial state to the target state in DFS. In other 
words, we can design A such that £ H DFS contains only 
the target state. We emphasis that although the con- 
trol field fjg{t) was specified to cancel Tr{C{p)A) in V, it 
makes contribution to the dynamics of the open system. 





FIG. 2: Fidelity of the system in the target state jDi) (a), 
and in the DFS (b). The control field fait) is turned off, i.e., 



/3(t)=0. n = 5,<^=f,/33 = f,7i 
Ao = 4, Ai = A2 = 2 and 7 = 1. 
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FIG. 1: The schematic energy diagram. A four-level system 
with two degenerate stable states |1) and |2) in external laser 
fields. The two degenerate states are coupled to the excited 
state |0) by two separate lasers with coupling constants f^i 
and ^2, respectively. While the stable state |3) is isolated 
from the other levels. The excited state |0) decays to \j) 
(j — 1, 2, 3) with decay rate 7^. 



III. EXAMPLE 

As an example of the Lyapunov control strategy, we 
discuss below a four-level system coupling to two external 
lasers, as shown in Fig. [TJ The Hamiltonian of such a 
system has the form 



i?o = ^^ij')oi + (E%|o)oi + ^-^-), 



(7) 



j=o 



j=i 



where Qj (j 



1,2) are coupling constants. Without 
loss of generality, in the following the coupling constants 
are parameterized as l^i = 17 cos 4> and fl2 — i^ sin cf) with 
fi = y/nl + ri|. The excited state |0) is not stable, it 
decays to the three stable states with rates 71, 72 and 
73, respectively. We assume this process is Markovian 
and can be described by the Liouvillian, 



3 _ ^ 



1 



-.P'^ ^7 



(8) 



j=i 



with (T = |0)(j| and a^ — {a- )t. It is not difficult to 
find that the two degenerate dark states 



\Di) = cos0|2) -sin0|l), 
1^2) = |3), 



(9) 



of the free Hamiltonian Hq form a DFS. Now we show 
how to control the system to a desired target state in the 
DFS. For this purpose, we choose the control Hamilto- 
nian He — X]i=i fji^)Hj with 



Hi = 



1111 
1111 
1111 
1111 

H2 = \Di){D2\ + \D2){Di\ 

H3 = \0){D2\ + \D2){0\. 




FIG. 3: Fidelity of the control with A = |-D2)(-D2| - |-Di>(I>i 
(i.e., the target state is |-Di)) as a function of the initial states. 
The parameters chosen are the same as Fig. [2l The control 
strength k,2 — 1,k,3 — 15 were specified for this plot. 



We shall use Eq. ((S]) to determine the control fields 
{fn{t)}, and choose 



I*) = sin/3icos/33|0)-Fcos/3icos/32|l) 
+ cos/3i sin/32|2) -t- sin/3i sin/^sjS) 



fll) 



(10) 



as initial states for the numerical simulation, where /3i , /32 
and /33 are allowed to change independently. We should 
emphasis that it is difficult to exhaust all possible ini- 
tial states in the simulation, because for a 4-dimensional 
system, there are 15 independent real parameters needed 
to describe a general state, even for pure states, 6 real 
independent parameters are required. The initial state 
written in Ea. dTTj) omits all (three) relative phases be- 
tween the states |0), |1), |2) and |3) in the superposition, 
and satisfies the normalization condition. fi(t) here is 
specified to cancel the contribution of Tr[£(p)74] to V, 
this means that /i(t) = -z^£f|Mg^, /2(i) and Ut) are 

determined by Eq.([5]). 

We have performed extensive numerical simulation 
with the initial states Eq. PT|) . Numerical results are pre- 
sented in Figs. [2]|4l The control field fsit) plays an im- 
portant role in this scheme as Fig. [5] shows. Fig. [5] tells 
us that without the control field fsit), the open system 
can be driven into the DFS (with A given below), but 
it can not be manipulated into a definite state in DFS. 
The physics behind is the following. With the given A 
(see below), /i(i) is always zero, so Hi plays no role in 
the control. The only control that enters the system is 
/2(t)i?2- From Eq.lO, we find that f2{t) takes zero pro- 
vided p = x\D2){D2\ + (1 - a;)lL»i}(Di|, (where x > 0), 
leading to the above observation. When the control field 
fait) is turned on. The four- level system can be con- 
trolled to a desired state in DFS by properly choosing A. 








FIG. 4: Fidelity of the control versus initial states. The target 
state is |D2> (or A = -ID2X-D2I + |Di)(Di|). K.3 = 15, the 
other parameters chosen are the same as in Fig[5]are chosen 
for this plot. 



deed further numerical simulations confirm this predic- 
tion that the control /13 can drive the system into \D2) 
with almost perfect fidelity 99%. The fidelity of the open 
system in the target state depends on the strength K3 — k 
of the control f^{t), the dependence is plotted in FiglS] 
With large k, the system would asymptotically converge 
to the target state as FiglS] shows. As expected, the con- 
trol fields f2{t) and fz{t) tend to zero when the open 
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FIG. 6: Control field /2(i) and fs,{t) as a function of time. 
= 5,0= f,^i = f,'^2 = f,^3 = f,K2 = 1, and«;3 = 15. 
A = \D2){D2\ — \D\){Di\. f\{t) is zero in this scheme. 



system converges to the target state, see Figl6] 



FIG. 5: Fidelity as a function of the control strength K3 = k. 

(a)i = \D2){D2\ - \D^){Dil (b)l = -|D2){D2| + |Di){Di|. 
<?!> = f,/3i = f,/32 = f,/33 = f, and K2 = 1. The other 
parameters chosen are the same as Fig. [21 



For example, when A = Ai = \D2){D2\ - |L'i)(£'i|, the 
system can be controlled into \Di) (see Fig[3]), whereas 
A = —Ai can drive the system into I-D2) (see FigH]). 
Based on the formalism in Sec. [Ill A = ^1 together 
with the controls could drive the system to the eigen- 
state of Ai with smallest eigenvalue (namely, \Di)), and 
to ID2) with A = —Ai. As figure 4 shows, however, 
the fidelity is not 1 for some initial states, for example 
/3i = 0. The reason is as follows. Though the choice of 
A = —Ai benefits the target state I-D2), since \D2) is 
the eigenstate of A with smallest eigenvalue, the control 
^3 = (|0)(^2| + h.c.) does not favor the control target 
[Z?2), because ^3 couples the states |0} and I-D2), and 
0) decays to the three stable states equally. This ob- 
servation suggests that /13 — (|0)(-Dil -|- h.c.) instead of 
H3 could help the control when the target is \D2). In- 



IV. CONCLUSION 

In summary, we have proposed a scheme to manip- 
ulate an open quantum system in the decoherence-free 
subspaces. This study was motivated by the fact that for 
Lyapunov control, it is usually difficult to optimally con- 
trol the system from an arbitrary initial state to a given 
target state, this is due to the LaSalle's invariant princi- 
ple. Our present study suggests that it is possible to drive 
a quantum system to a desired state in DFS by elabo- 
rately designing the controls. The results do not break 
the LaSalle's role, instead it reduces the invariant set £ 
to include the target state only. To demonstrate the pro- 
posal we exemplify a four-level system and numerically 
simulate the controlled dynamics. The dependence of the 
fidelity on initial states as well as the control fields are 
calculated and discussed. This scheme put the Lyapunov 
control on quantum open system one step forward, and 
shed light on the quantum control in DFS. 
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